arXiv: 1505.05150v2 [cond-mat.quant-gas] 19 Oct 2015 


Dynamical Quasicondensation of Hard-Core Bosons at Finite Momenta 


L. Vidmar/’^ J. P. Ronzheimer,^’^ M. Schreiber,^’^ S. Braun,S. S. Hodgman,^’^ 

S. Langer,^ F. Heidrich-Meisner,^’^ I. Bloch,and U. Schneider^’® 

^Arnold Sommerfeld Center for Theoretical Physics, 
Ludwig-Maximilians-Universitdt Munchen, D-80333 Munchen, Germany 
^Department of Physies, Ludwig-Maximilians-Universitdt Munchen, D-80333 Munchen, Germany 
^Max-Planek-Institut fiir Quantenoptik, D-85748 Garching, Germany 
'^Department of Physies and Astronomy, University of Pittsburgh, Pittsburgh, Pennsylvania 15213, USA 

® Cavendish Laboratory, University of Cambridge, 

J. J. Thomson Avenue, Cambridge CB3 OHE, United Kingdom 

Long-range order in qnantum many-body systems is usually associated with equilibrium situa¬ 
tions. Here, we experimentally investigate the quasicondensation of strongly-interacting bosons at 
finite momenta in a far-from-equilibrium case. We prepare an inhomogeneous initial state consist¬ 
ing of one-dimensional Mott insulators in the center of otherwise empty one-dimensional chains in 
an optical lattice with a lattice constant d. After suddenly quenching the trapping potential to 
zero, we observe the onset of coherence in spontaneously forming quasicondensates in the lattice. 
Remarkably, the emerging phase order differs from the ground-state order and is characterized by 
peaks at finite momenta ±(7r/2)(li/d) in the momentum distribution function. 
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The nonequilibrium dynamics of quantum many-body 
systems constitutes one of the most challenging and in¬ 
triguing topics in modern physics. Generically, interact¬ 
ing many-body systems are expected to relax towards 
equilibrium and eventually thermalize [1, 2]. This stan¬ 
dard picture, however, does not always apply. 

In open or driven systems, one fascinating counterex¬ 
ample is the emergence of novel steady states with far- 
from-equilibrium long-range order, i.e., order that is ab¬ 
sent in the equilibrium phase diagram. This includes 
lasers [3] , where strong incoherent pumping gives rise to a 
coherent emission, and driven ultracold atom systems [4] . 
The emergence of order far from equilibrium is also stud¬ 
ied in condensed matter systems [5] and optomechanical 
systems [6]. 

In recent years, closed quantum systems without any 
coupling to an environment have come into the focus of 
experimental and theoretical research. Experimental ex¬ 
amples range from ultracold atoms [7-15] to quark-gluon 
plasmas in heavy ion collisions [16]. In closed, non-driven 
systems, two famous examples for the absence of thermal- 
ization [8, 9, 14] are many-body-localized [17-19] and in- 
tegrable systems [20]. These peculiar systems allow for 
nonergodic dynamics and novel quantum phenomena. 

Spontaneously emerging order is in general associated 
with equilibrium states at low temperatures. The canoni¬ 
cal example is the emergence of (quasi-) long-range phase 
coherence when cooling an ideal Bose gas into a Bose- 
Einstein (quasi-) condensate [21, 22]. In this case, ther¬ 
modynamics ensures that, for positive temperatures [23], 
the single-particle ground state becomes macroscopically 
occupied and thereby dictates the emerging order. Even 
in studies of the nonequilibrium dynamics at quantum 
phase transitions [24] , the emergence of coherence is typ¬ 
ically associated with gently crossing the transition from 
an unordered into an ordered state, and the strongest 



FIG. 1. Quasicondensation of bosons, (a) Density of states 
p(e) of a homogeneous ID lattice, (b) Dispersion e{q) (solid 
line) and group velocity Vg{q) (dotted line) versus quasimo¬ 
mentum q. (c) Sketch of quasimomentum distribution n{q): 
In equilibrium, ID bosons quasicondense at the minimum of 
the band at q = 0, while in a sudden expansion the quasicon¬ 
densation of hard-core bosons occurs at hq = ±{n/2){h/d). 
This quasimomentum lies in the middle of the spectrum and 
is consistent with the vanishing energy per particle of this 
closed many-body system. 


correlations and largest coherence lengths appear in the 
adiabatic limit [25, 26]. 

Here, in contrast, we study a condensation phe¬ 
nomenon of strongly interacting lattice bosons far from 
equilibrium. After a sudden quantum quench, we exper¬ 
imentally observe the spontaneous emergence of a long- 
lived phase order that is markedly different from the 
equilibrium order (cf. Fig 1). To this end, we prepare a 
density-one Mott insulator of strongly interacting bosons 
in the center of a three-dimensional (3D) optical lattice. 
Next, we transform the system into an array of indepen¬ 
dent one-dimensional (ID) chains, entering the regime 
of integrable hard-core bosons (HCBs). By suddenly 
quenching the conhning potential along the chains to 
zero, we induce a sudden expansion of the cloud in a ho- 
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mogeneous lattice [12, 13, 15, 27] with a lattice constant 
d and detect the formation of a non-ground-state phase 
profile as a dynamical emergence of peaks at momenta 
hk = zL{7r/2){h/d), half way between the middle and 
the edge of the Brillouin zone, in time-of-flight (TOF) 
distributions. This finite-momentum quasicondensation 
was first discussed by Rigol and Muramatsu [28] (see also 
Refs. [29-33]), but has not been studied experimentally 
so far. 

Ideal case. The idealized setup to study finite- 
momentum quasicondensates is shown in Fig. 2. We 
consider the Hamiltonian H = —J + h.c.), 

where creates a HCB on site j of a ID lattice. The 
infinitely large on-site repulsion is accounted for by the 
hard-core constraint (d])^ = 0. The initial state is a 

product state \ijjo) = HjeLo completely filling the 

central Lq sites of an otherwise empty and infinitely large 
ID lattice. This initial state consists of iV = Lq localized 
particles with a flat quasimomentum distribution and 
contains no off-diagonal correlations, i.e., (d|dj+r) = 0 
for r 7 ^ 0. Surprisingly, the quasimomentum distribution 
n{q) = j- J2j I develops singularities at fi¬ 

nite quasimomenta hq = ±(7r/2)(d/d) during the expan¬ 
sion (^E > 0). As shown in Fig. 2(c), these singularities 
correspond to the emergence of power-law correlations 

(dldj+r) = A(r)e*‘^(’'); A{r) ~ r"5; $(r) = ±^r (1) 

in each half of the expanding cloud, shown in Fig. 2(d). 
These power-law correlations justify the name quasicon¬ 
densate [28]. Curiously, the exponent 1/2 equals the 
ground-state exponent [28, 30], even though the system 
is far away from equilibrium, with the energy per parti¬ 
cle being much higher than in the ground state. In con¬ 
trast to the ground state, the correlations show a running 
phase pattern $(r) with a phase difference of ±7r/2 be¬ 
tween neighboring lattice sites, giving rise to peaks at 
finite quasimomenta. Coherence and quasicondensation 
emerge independently in the left- and right-moving halves 
of the cloud, corresponding to two macroscopically occu¬ 
pied degenerate eigenstates of the one-particle density 
matrix (d^d;) that have spatial support in the left- or 
right-moving cloud, respectively [28]. This quasiconden¬ 
sation at finite quasimomenta can equivalently be seen as 
quasicondensation at q = 0 in the respective co-moving 
frames [32]. 

In one dimension, HCBs can be exactly mapped to 
noninteracting spinless fermions via the Jordan-Wigner 
transformation [34] . By virtue of this mapping, the den¬ 
sity nj = {a^jdj) of HCBs is identical to that of free 
fermions for all times, whereas the same is not true for the 
quasimomentum distribution [35, 36]. While the occupa¬ 
tions of fermionic quasimomenta are constants of motion, 
the non-local phase factors in the Jordan-Wigner trans¬ 
formation give rise to the intricate momentum dynamics 
studied here. In the ideal case, the dynamical quasicon¬ 
densates form over a time scale 1^ ^ Q.3Nt [28, 33], 
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FIG. 2. Finite-momentum quasicondensation under ideal 
conditions [28]. We consider A = 50 initially localized HCBs. 
(a) Density Uj as a function of time, (b) Density as a func¬ 
tion of the rescaled coordinate x = {j — 25.5) / {2tE /t) in the 
region bounded by the dashed lines in (a). For i < 1, the 
data collapse to the scaling solution [39] n(i) = arccos {x)/tt. 

(c) Quasimomentum distribution n{q) as a function of time. 

(d) One-particle correlations at ts = 0.24Ar. Main panel: 
|(a]aj+r)| at j = 26 (circles) and A{r) = aj^/r (line) with 
a = 0.29 [40]. Inset: Phase pattern <I>(r). 


where N is the number of particles in the initial state 
and T = h/ J denotes the tunneling time. For very long 
times, n(g) slowly decays back to its original flat form 
as a consequence of the dynamical fermionization mech¬ 
anism [33, 37, 38]. 

The dynamical quasicondensation at finite momenta 
is an example of a more general emergence of coher¬ 
ence in a sudden expansion. For instance, interacting 
fermions described by the Fermi-Hubbard model exhibit 
ground-state correlations in the transient dynamics as 
well [32]. Furthermore, there is a close connection to 
quantum magnetism, as can be seen by mapping HCBs 
to a spin-1/2 XX chain: The transient dynamics in each 
half of the expanding cloud of HCBs is equivalent to the 
melting of a domain-wall state [39, 41-49] of the form 
[■i/o) = It-•-ttii- ■ -i)- For this problem, a scaling so¬ 
lution exists [39], which also applies to the sudden ex¬ 
pansion at times t < As a consequence, the densi¬ 
ties Uj measured at different times collapse onto a single 
curve, as shown in Fig. 2(b). Furthermore, for spin-1/2 
XX chains the emergence of power-law decaying trans¬ 
verse spin correlations modulated with a phase of 7r/2 
has been derived analytically [44]. An interesting per¬ 
spective onto the emergence of coherence results from 
noticing that both our expanding bosons and the melt¬ 
ing domain-walls realize current-carrying nonequilibrium 
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FIG. 3. Dynamics during sudden expansion and time-of-flight sequence, (a) Sketch of the experimental sequence: We start 
with a trapped gas that is released from the initial confinement at t = 0 (top) and then expands for a time ts in the optical 
lattice (middle). At Ie, all potentials are removed and the distribution is measnred after a hnite time-of-flight time txoF 
(bottom), (b)-(e): In-situ density distributions during the sudden expansion, integrated along the y- and a-axes. The lattice 
constant corresponds to d = A/2 « 368nm. (f)-(i): TOF density distributions taken at txoF = 6ms for the same Ib as in 
(b)-(e). All measured density distributions are averaged over 9 to 11 experimental realizations. 


steady states (see [50] for a discussion). 

Experimental set-up and results. The experimental set¬ 
up is identical to that employed in our previous experi¬ 
ment on in-situ density dynamics [13]. We load a Bose- 
Einstein condensate of approximately 10® atoms 
from a crossed optical dipole trap into a blue-detuned 
3D optical lattice with a lattice depth of Vb « 20Er, 
where E^ = h^/{2mX^) denotes the recoil energy with 
atomic mass m and lattice laser wavelength A « 737 nm. 
During the loading of the lattice, we use a magnetic Fes- 
hbach resonance to induce strong repulsive interactions 
between the atoms, suppressing the formation of doubly 
occupied sites [13]. This results in a large density-one 
Mott insulator in the center of the cloud. We hold the 
atoms in the deep lattice for a 20 ms dephasing period, 
during which residual correlations between lattice sites 
are mostly lost such that the atoms essentially become 
localized to individual lattice sites [51]. The expansion is 
initiated at ts = 0 by simultaneously lowering the lattice 
depth along the expansion axis in 150 /iS to Vq ^ 8 E^ 
(setting J Ki h X 300Hz, r « 0.5 ms), and reducing the 
strength of the optical dipole trap to exactly compensate 
the anti-confinement of the blue-detuned lattice beams. 
This creates a flat potential along the expansion direc¬ 
tion. Figures 3(b)-3(e) show the ballistic expansion of 
the in-situ density [13], monitored using absorption imag¬ 
ing. During the deep lattice period, the magnetic field 
is changed to tune the on-site interaction strength dur¬ 
ing the expansion to U = 20J. We have numerically 
verified that the essential features of dynamical quasi¬ 
condensation are still present for this value of U / J, with 
a shift of the peak position by « 10% towards smaller 


values [31, 50]. 

In order to measure the momentum distribution as a 
function of expansion time t-^, we employ an adapted 
TOF imaging technique. Directly before shutting off 
all lattice and trapping potentials, we rapidly increase 
the lattice intensity along the expansion axis for 5 
to a depth of 3375^. This time is too short to affect 
correlations between different sites and the momentum 
distribution. Nonetheless, it results in a narrowing of 
the Wannier functions, which leads to a broadening of 
the Wannier envelope in the TOF density distribution 
and thereby facilitates the observation of higher-order 
peaks [50]. 

Figures 3(f)-3(i) contain the main result of our experi¬ 
ment, namely the TOF density distributions, which cor¬ 
respond approximately to the momentum distribution, 
taken at different expansion times Ie- In Fig. 3(f) the 
TOF sequence was initiated at Ie = 0, i.e., directly after 
initiating the expansion. We observe a central peak at 
k = 0 and two higher order peaks at Hk = ±{2TT){li/d), 
indicating a weak residual k = 0 coherence probably re¬ 
sulting from an imperfect state preparation. During the 
expansion, however, the momentum distribution changes 
fundamentally and the remnants of the initial coherence 
quickly vanish. Instead, new peaks at finite momenta 
are formed. These peaks directly signal the sponta¬ 
neous formation of a different phase order. This is best 
seen in Fig. 3(i) at Ie = 36t, where the finite-momenta 
peaks are clearly established. The observed peak posi¬ 
tions correspond to the expected momenta close to hk = 
±{Tr/2){h/d) [50]. In addition, Figs. 3(g) and 3(h), taken 
at tE = 7r and I4r, respectively, hint at a hne structure 
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of the emerging peaks. This structure is a consequence 
of the finite TOF time tTOF = 6ms, which results in the 
TOF distributions being a convolution of real-space and 
momentum-space densities. We sketch this situation in 
Fig. 3(a). As discussed before, the peaks in n{k) close to 
hk = —{'K/2){h/d) and -\-{'k/ 2){h/d) originate from the 
left- and right-moving portions of the cloud, respectively. 
Due to the finite tTOF, the higher-order peak of the left- 
moving cloud with momentum (—7r/2 -|- 2tt) (h/d) and 
the main peak of the right-moving cloud with momen¬ 
tum {tt/ 2){h/d) (and vice versa) may overlap in the TOF 
data. A perfect overlap gives rise to single sharp peaks 
such as the ones present in the data shown in Fig. 3(i), 
while shorter expansion times, as shown in Figs. 3(g)- 
3(h), result in a partial overlap and additional structure 
(see [50] for details). 

Comparison with exact time evolution. We numerically 
model the dynamics of ID HCBs for realistic conditions: 
(i) The experimental set-up consists of many isolated ID 
chains, which are not equivalent due to the 3D harmonic 
confinement. Experimentally we can only measure an en¬ 
semble average over all tubes, (ii) Both the finite temper¬ 
ature of the original 3D Bose-Einstein condensate as well 
as nonadiabaticities during the lattice loading result in a 
finite entropy, and thereby holes, in the initial state. We 
therefore average the results over different initial prod¬ 
uct states drawn from a thermal ensemble of a harmon¬ 
ically trapped 3D gas of HCBs in the atomic limit [50]. 
Chemical potential and temperature were calibrated to 
reproduce the experimental atom number and an average 
entropy per particle of 1.2 fcs [52], thereby leaving no free 
parameters for the simulations. To test the consistency of 
the approach, we compare the average density nj during 
the expansion with the in-situ images in Figs. 3(b)-3(e) 
and find a good agreement. In addition, the time evo¬ 
lution of the half width at half maximum of the density 
distribution, shown in Fig. S4 of [50], is consistent with 
the ballistic dynamics as previously measured in the same 
experimental set-up [13]. 

Since the momentum distribution is experimentally 
measured at a finite ItoF) we explicitly calculate the 
TOF density distributions without employing the far- 
held approximation [50, 53] and compare the results to 
the experimental data in Figs. 3(f)-3(i). Remarkably, 
the positions and the structure of the peaks agree very 
well between experiment and theory, thereby supporting 
our two main results: The central peaks indeed corre¬ 
spond to a large occupation of quasimomenta close to 
hq — ±{-K/2){h/d), i.e., to a bunching of particles around 
the fastest group velocities in the middle of the single¬ 
particle spectrum. In addition, the hne structure visi¬ 


ble for intermediate expansion times [cf. Figs. 3(g)-3(h)], 
which becomes more apparent in the experiment when 
comparing different Itof (cf. Fig. S2 in [50]), directly 
conhrms the independent emergence of coherence in the 
left- and right-moving portions of the cloud. Compared 
to the ideal case, the presence of holes in the initial state 
causes a reduced visibility of the TOF density distribu¬ 
tions [50]. Moreover, the finite initial entropy gives rise 
to a crossover of one-particle correlations from a power- 
law decay at short distances to a more rapid decay at 
long distances [50], similar to the effect of a finite tem¬ 
perature [54] in equilibrium. 

We attribute the discrepancies between experimental 
and numerical results at short times, see Fig. 3(f), to the 
weak residual fc = 0 coherence in the initial state. Addi¬ 
tional discrepancies may arise because of a small admix¬ 
ture (< 5%) of doublons in the initial state [13, 15] as well 
as small residual potentials, yet we conclude that these 
play a minor role [50] . Compared to the previously stud¬ 
ied time dependence of density distributions and expan¬ 
sion velocities [13], the momentum distribution is more 
sensitive to such imperfections [50] . Performing a similar 
experiment with a single ID system would allow the pre¬ 
dicted scaling of and the maximum peak height with 
atom number [28] to be experimentally tested. 

Conclusions and outlook. We have reported experi¬ 
mental evidence for a far-from-equilibrium quasiconden¬ 
sation at finite momenta of expanding ID HCBs in an 
optical lattice. The expanding particles bunch at quasi¬ 
momenta close to zt{Tr/2){h/d) and the analysis of TOF 
distributions demonstrates the existence of two indepen¬ 
dent sources of coherence. 

Whether such dynamical condensation persists in 
higher dimensional systems constitutes an open prob¬ 
lem, given that the existing theoretical results are based 
on exact diagonalization of small systems [55] or time- 
dependent Gutzwiller simulations [56]. Both future ex¬ 
periments or advanced numerical methods (see, e.g., [57, 
58]) could help clarify this question. More generally, our 
results raise the question of whether this type of sponta¬ 
neously emerging coherence is limited to integrable sys¬ 
tems and whether genuinely far-from-equilibrium order 
can also occur in generic closed many-body systems. 
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FIG. SI. Experimental TOE density distributions. Compar¬ 
ison between a measurement (a) without additional lattice 
pulse and (b) with additional pulse before the TOF sequence 
(in both cases Ie = 22r and Itof = 12ms). Vertical lines are 
guides to the eye. Density distributions are averaged over 10 
experimental realizations. 


Supplemental material 


SI. EXPERIMENTAL DETAILS 
A. Increasing the visibility of higher order peaks 

By switching off the lattice instantaneously at the be¬ 
ginning of the time-of-flight (TOF) sequence, each quasi¬ 
momentum Hq is projected onto a superposition of states 
with free-space momentum hk (with k = q) and higher- 
order (Bragg) contributions with momenta fi,(fc±n27r/d), 
n € N. In order to increase the visibility of the higher- 
order peaks, we increase the intensity of the cc-lattice for 
5 pLS to approximately 33i?r immediately before releas¬ 
ing the atoms from the lattice. This time is long enough 
to lead to a narrowing of the on-site (Wannier) wave- 
functions, resulting in a broadening of their momentum 
distributions, while being short enough to have no signifi¬ 
cant influence on the quasimomentum distribution. This 
leads to a broadening of the envelope of the momentum 
distribution observed after time of flight and thereby en¬ 
hances the visibility of the higher order (Bragg) peaks, 
as shown in Fig. SI. 


B. Analysis of observed peak positions during time 
of flight 

The expected signatures of the quasicondensation at 
finite quasimomenta close to ±(7r/2)(Ii/(i) are distinct 
peaks in the density distribution after the TOF sequence. 
Since the quasicondensates form independently in the 
left- and right-moving portions of the cloud during the 


sudden expansion, the peaks observed in the TOF distri¬ 
bution can be interpreted as the sum of two independent 
contributions originating from different points in space 
denoted by d!^ and respectively. For specific TOF 
times, different momenta stemming from different posi¬ 
tions in the cloud can become superimposed and appear 
as one larger peak in the measured TOF density distri¬ 
bution. To demonstrate this effect, we extract the peak 
positions for various different TOF times ixOF with a 
simple fit function consisting of six Gaussian peaks on 
top of a broad background. 

The six small Gaussian peaks are grouped into two 
groups corresponding to the left- and right-moving por¬ 
tions of the cloud. Each group contains one peak at po¬ 
sition cZl/r that corresponds to d^^^ + t^OF ■ 
respectively, and the additional higher order peaks at po¬ 
sitions dL/R±d 27 r, where fi 27 r = Itof ‘^Trh/(dm). In order 
to further reduce the number of free fitting parameters, 
we assume all six Gaussian peaks to have the same width 
wc, but allow for individual amplitudes ^r/r- com¬ 
plete fit function F(x ,...) is given by 

F{x ,...) = Pbix, Ab, db, Wb) 

+ Pf{Xt , A^^ A\^^di^,wa) (SI) 

+ dR, Wg) 

with the broad background peak 

Pb{x,...)= Abe\ 


and the groups of small peaks 
PF/Ki.x, ...)=A ^ 


■AO 


L/R 

l/rC 
1 

L/R 


(®-^L/R“^27r)^ 

e V 


(z-dL/R)- 


(z-dL/R + d27r)^ 
-2- 


(S3) 


■ AJ /dC ^ '“"G 


Due to the additional sudden switching off of the mag¬ 
netic field used to address the Feshbach resonance and 
the resulting eddy currents, there are additional position- 
dependent forces acting on the atoms during the initial 
stage of the TOF sequence. In order to correctly map 
the observed distances after time of flight onto the initial 
momenta, we apply a correction that is determined by 
monitoring the evolution of the TOF peaks of an equilib¬ 
rium q = 0 Bose-Einstein condensate released from the 
lattice, where the momentum composition is known to 
consist of peaks at ±n2TT{h/d) with n S N. 

Figures S2(a)-S2(d) show the density distributions af¬ 
ter an expansion time ts = 36r in the lattice and varying 
TOF times ixoF together with the resulting fits. Even 
though the fits result in strongly varying peak ampli¬ 
tudes, they nonetheless serve to faithfully identify the 
peak positions. 

In Fig. S2(e), we show the extracted positions of the 
six peaks (red and blue circles) for a fixed expansion time 
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tjoF = 6ms tjof = 10ms tjof = 12ms tjop = 15ms 




FIG. S2. Peak positions as a function of time-of-flight time. 
(a)-(d) Density distributions obtained for a fixed expansion 
time of Ie = 36r (blue solid lines). Black dashed lines are fits 
with a sum of six Gaussian peaks and a broad background, 
black solid lines show only the contributions of the six peaks. 
All measured density distributions are averaged over 9 to 10 
experimental realizations, (e) Extracted peak positions for 
Ie = 36t (dots). Dashed lines are linear fits, solid black lines 
show the expected slope for peaks traveling with momentum 
±{'n/2){h/d). Error bars in the position are smaller than the 
data points. 


Ie = 36t. The red circles are associated with the compo¬ 
nent at negative momentum hki^ and its corresponding 
higher order peaks at h{ki, ± l-nld)., the blue circles are 
associated with the positive momentum hk^ and higher 
order peaks. 

The dashed lines are linear fits, their slopes represent 
the extracted velocities. For comparison, the solid black 
lines show the theoretically expected slope for momenta 
±{Tr/2)(h/d). The peaks do indeed propagate with mo¬ 
menta close to the ideal expectation; the small shift to¬ 
wards smaller momenta can be attributed to the finite 
interaction strength U/J < oo used in the experiment 
(see Sec. S3 C). 

The extrapolation to Itof = 0 shows that the two 
momentum groups indeed originate from two different 
points in space in the initial cloud, with the left-moving 
part starting further to the left and vice versa. 


S2. TIME-OF-FLIGHT DISTRIBUTION 
WITHOUT THE FAR-FIELD APPROXIMATION 


After the lattice flashing, each site of the lattice can 
be described by the Wannier function wo{r), which we 
approximate by a Gaussian of width a: 

wo{r) = —(S4) 

y/cry/TT 

As a consequence of the lattice flashing method de¬ 
scribed in Sec. SI A, a ^ d. After switching off the 
lattice, each Wannier function freely evolves in time. We 
calculate the free time evolution in momentum space, 
FT[ryo(r,tTOp)] = FT[r(;o(r)]e“*^''‘'^°^/^, where FT de¬ 
notes the Fourier transform and Ek = (hkY/(2m) is the 
free-space single-particle dispersion. For the propagation 
in real space, we get 


Wo(r,tTOF) = 




• Mtof 


)V^ 


(S5) 


The total density profile nTOF(F) = ?t-tof(g ^tof) = 
{'ijd (r)'i/i{r)) at a given fxOF is a sum over contributions 
from all (occupied) lattice sites, 

riTOFir) = ^ Wo(r - , <tof)wo(?' - n,iTOF)(4.dn), 

Tj.ri 

(S6) 

where the field operator ip(r) was expanded in time- 
evolved Wannier functions. In the case of one¬ 
dimensional (ID) lattice dynamics studied here, off- 
diagonal correlators (a]..dr/) are nonzero along a single 
direction only, hence we define x = r/d and /j, = 

The density can be written in a compact form by intro¬ 
ducing two dimensionless parameters a = l/[(a/d)'^ -\- 
(ht'TOF/(‘Fnad))'^] and (3 = ahtTOF/("mcr^), resulting in 

nTOF(x) = x 

^ j.l 

This generalizes the result of Gerbier et al. [53], where 
only the phase factors have been taken into account. In 
our experiment cr = 50 nm, i.e., cr/d « 0.13. The central 
task is therefore to calculate the correlators (ajfi,;), for 
which we provide details in the following Secs. S3 A and 
S3B. 


S3. DETAILS OF THE NUMERICAL 
SIMULATIONS 


At Itof = 6ms in Fig. S2(e), we can also observe how 
multiple peaks can come to overlap during the TOF evo¬ 
lution, resulting in the observation of fewer but stronger 
peaks, explaining the number of peaks observed in, e.g.. 
Fig. S2(a). 


A. Realistic modeling of the initial state 

We model the experimental set-up by calculating the 
exact time evolution of many ID systems of hard-core 
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FIG. S3. Distribution of probabilities P{N) to find par¬ 
ticles in chains with particle number N. We compare the 
zero-entropy {S = 0) and the finite-entropy case {S/N 30 = 
1.2 fcfl). We use a binning with AN — 10. 



FIG. S4. Time dependence of the core radius Rc{tE) at 
S/N 30 = 1.2 fcs. The solid line represents the fit function 
Rcits) = UctB+7 with Vc = 2.06(d/T) and 7 = 14.8, obtained 
by fitting to the numerical data in the time interval t^/r € 
[30,50]. 


bosons (HCBs) [28, 33] that differ by particle numbers 
N and their distribution in the initial state. The generic 
initial state within a single chain is a product state 

iv^o)=n 

jeJ 

where ff does not necessarily represent a sequence of 
neighboring sites. The initial states are drawn from a 
thermal distribution of HCBs in the atomic limit (J = 0) 
of a 3D cloud and the final result of the simulations is 
obtained by summing up the data for initial states with 
different J'. 

The probability to find a particle at a given lattice site 
is given by the following partition function, 

Z= ^ ^-P{nV{x,y,z)-nno) ^ 
n=0,l 

where fj-o represents the chemical potential, j3 = l/ifk^eT) 
is the inverse temperature and the harmonic confinement 
is given by 


V{x,y,z) = 


mcuz 


[x^ + y^ + iAzf], 


(SIO) 


with 0Jx,y = w = 27r X 60 Hz and the aspect ratio 
A = ojzjujx^y = 2.63. The ID systems are aligned along 
the x-direction and the results are integrated over chains 
at different positions (y, z). Due to the harmonic confine¬ 
ment, chains at different transverse positions can have 
vastly different atom numbers in the range N G [0,80]. 
Clearly, for T —i 0, the initial state in the 3D cloud con¬ 
sists of an ellipsoid of singly-occupied sites without any 
holes. A finite-entropy initial state, on the other hand, 
contains a finite density of holes that increases when mov¬ 
ing outwards. Equation (S9) leads to the global density 


n(x,y,z) 


1 

gp(V(x,y,z)-uo) -(- 1 ■ 


(Sll) 


We determine f and /tq by matching the total particle 
number A^ 3 d of the 3D cloud and the average entropy per 


particle S/N^u to experimental conditions. The average 
entropy per particle is defined as 

S/N 3 u= Y. SJN 31 ,, (S12) 

i={x,y,z) 


where 

Si = -ks (uilogn* + (1 - n01og(l - n*)) (S13) 

and rii represents the probability to find a particle at site 
i. In our experiment, we have A 3 D = y z '>T-ix, y, z) « 
0.9 X 10^ and we estimate S/N^d « 1.2 ks [52]. A typi¬ 
cal configuration for = 50 is shown in Fig. S5(a2). 
Note that we also symmetrize each configuration by tak¬ 
ing its mirror image (with respect to the central site of 
the originally occupied region). 

In Fig. S3 we plot the distribution of probabilities 
P{N) to find particles in chains with a given number 
of particles N. This information reveals the relative con¬ 
tribution of chains with a given particle number to the 
dynamics. The distribution for S' = 0 is shown with cir¬ 
cles. In this case, the 3D state corresponds to a perfectly 
filled ellipsoid Rq = x^ y"^ (Xz)^ with Ro/d = 38 

and, as a consequence, P{N) oc for N < 70. As a 
general trend, the distribution shifts towards smaller N 
with increasing entropy (squares in Fig. S3). 

In our calculations, we simulate the dynamics of 
representative chains with particle numbers N G 
{1,10,20,30,40,50,60,70} and the relative weights 
shown in Fig. S3 for «S'/A^ 3 d = 1.2 ks- All together we 
average the data over 100 different initial configurations 
and we have verified that the sampling over this subset 
of all possible N is sufficient. In all data for in-situ and 
TOF density distributions shown in the figures, we con¬ 
volve the numerical data with a Gaussian filter of width 
Id to account for the finite experimental imaging resolu¬ 
tion. 

Figure S4 shows that the core radius (i.e., the half 
width at half maximum) increases linearly in time ac¬ 
cording to Rc{t-E) = VctE + const with the core velocity 
very close to Vc = 2{d/T), in agreement with [13]. This 
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(a) Initial state 


(a1) 

(a2) 


. ... I I I . I I I... I .. . I ■ > I I I I.. . I 

■V..VV.VVV..VV..VVV.VVV.VV.V.■VVVV.^VVVVVVVVVVV..VVVVV.VVVVVV.■S'.VVV..V. . .V. 


(b) Time evolution 





FIG. S5. Influence of holes in the initial state on the one-particle correlations, (a) Snapshots of atom configurations in the 
initial state {N = 50 in all cases) for (al) an ideal initial state I'i/'o), Eq. (S14), and for (a2) a typical initial state in Eq. (S8), 
generated at S/A^sd = 1-2 ks- (b) Modulus of one-particle density matrix |(d^ai)| for (bl) Ie = lOr, emerging from the ideal 
initial state in (al), and for (b2) t-E = 14r, emerging from the symmetrized initial state consisting of the state in (a2) and 
its mirror image, obtained by reflection with respect to the central site. (b3) Spatial decay of |(o]aj+r)| at j = 68 of the 
(non-symmetrized) state shown in (a2) at t^ = 14t. The solid line is the function /(r) = (with a = 6.1 x 10“®) to 

illustrate that the decay at large distances is faster than a power law. 


suggests that the presence of holes in the initial state 
does not affect the ballistic expansion of ID HCBs. 


B. Correlation functions from the exact time 
evolution of a 3D cloud of HCBs 

We proceed with the analysis of the one-particle corre¬ 
lations (a]dz) during the time evolution. We compare a 
representative chain at S/N^-q = 1.2 ks to the ideal ini¬ 
tial state, both with TV = 50 particles. The ideal initial 
state is (see the main text) 

iV'o) = n 

leio 

where the length of the region Lq equals N. The distri¬ 
bution of particles in \tpQ) is displayed in Fig. S5(al), 
while the particle configuration of a typical chain at 
S/N^j:) = 1.2 ks is shown in Fig. S5(a2). 

The modulus of the one-particle density matrix | (a^a;) | 
resulting from |'0o) is shown in Fig. S5(bl) for Te = lOr. 
The two squares visualize the two sources of coherence 
emerging in the left- and right-moving cloud [28], corre¬ 
sponding to the peaks in n{q) at hq = —{TT/2){H/d) and 
hq = {-K/2){h/d) in Fig. 2(c), respectively. Within the 
squares, the correlations decay according to a power law 
(cf. Fig. 2(d) in the main text). 


The modulus of the one-particle density matrix j (a^d/) j 
for the chain displayed in Fig. S5(a2) is shown in 
Fig. S5(b2) at Te = 14t. As this example illustrates, 
the correlations in chains with holes show a less regular 
pattern, compared to the ideal situation. However, the 
calculated phase differences between neighboring sites ex¬ 
actly equal ±7r/2, identical to the ideal case shown in the 
inset of Fig. 2(d) of the Letter. In Fig. S5(b3) we plot 
|(d^dj+r)| as a function of r at j = 68 and the compar¬ 
ison to the function f(r) = clearly shows that, 

taking holes in the initial state into account, the cor¬ 
relations decay much faster than in the ideal case for 
large distances. Therefore, even though signatures of the 
finite-momentum quasicondensation can still be detected 
in a system with a finite initial entropy, the correlations 
have become short ranged. This is equivalent to the be¬ 
havior of equilibrium ID Bose gases at any finite temper¬ 
ature [59], where temperature cuts off the ground-state 
power-law decay. The temperature, or the amount of en¬ 
tropy, thereby sets the crossover scale from power-law to 
exponentially decaying correlations. 

C. Other sources for a loss of coherence 

While the overall agreement between the experimen¬ 
tal data and our numerical simulations is quite good 
(see Figs. 3(f)-3(i) of the Letter), there are neverthe- 
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FIG. S6. Quasimomentum distribution of HCBs and in¬ 
teracting bosons at U/ J — 20. We compare the dynamics 
at Ie = {N/4:)t and N = 16. Results for HCBs are shown 
with dashed lines in both panels. In (a), we show n{q) for 
the expansion from the ground state (g.s.) at U/J — 20 and 
the initial product state (p.s.) \ipo) from Eq. (S14). The solid 
and dashed vertical lines denote the estimate q' = 1.42/d 
and 1.47/d, respectively, for the quasimomenta at the shifted 
maxima. In (b) we show n{q) for the expansion from an ini¬ 
tial product state with a single doublon and hole (dash-dotted 
line; averaged over 184 configurations) and with two doublons 
and two holes (solid line; averaged over 552 configurations). 


less small discrepancies. For instance, the experimental 
data exhibits broader peaks and a lower visibility even 
at the longest expansion times [see Fig. 3(g) in the main 
text]), which arise from additional sources of decoherence 
present in the experiment. We have numerically investi¬ 
gated the effect of the following deviations from ideal con¬ 
ditions: (i) the finite interaction strength U/J = 20<oo 
used in the experiment, (ii) the presence of small resid¬ 
ual potentials, (iii) the presence of a small admixture of 
doublons in the initial state, and (iv) the influence of 
coherence in the initial state, i.e., the presence of a max¬ 
imum in 'n{q) at q = 0. Away from the HCB limit [i.e., 
investigating the effects (i) and (iii)], one needs to resort 
to time-dependent density matrix renormalization group 
(tDMRG) simulations [60-62]. For the tDMRG method, 
a realistic modeling with particle numbers as large as in 
the experiment and sufficiently long times is currently 
unfeasible. 

The dynamics of the quasimomentum distribution in 
the sudden expansion at U/ J < oo starting from the cor¬ 
related ground state at density one has been investigated 
in [31]. It has been observed that the position of the 
maxima moves to smaller values of li\q\ < [tt/ 2){h/d). 
As a heuristic estimate for the peak position [28] and to 


account for the conservation of total energy in this quan¬ 
tum quench, one can assume that all particles condense 
at the point in the single-particle dispersion correspond¬ 
ing to the average energy per particle. This results in the 
condition E/N = —2Jcos{qd). In the hard-core case, 
the average energy per particle vanishes, and hence one 
obtains hq = ±{Tr/2)(h/d). 

For expansion from the correlated ground state at 
U/J = 20, we express the ground-state energy Ag.s, = 
—ANJ'^/U by second-order perturbation theory around 
the J/U = 0 limit, arriving at hq' = aiccos {2J/U){h/d). 
Our tDMRG results for n{q) at U/J = 20 are shown in 
Fig. S6(a) and the estimate for the position of the peak 
hq' = arccos (2/20)(/i/d) = lA7{h/d) is in good agree¬ 
ment with the numerical data [31]. 

In our experiment, we ideally start from the product 
state (S14) of one boson per site. The quasimomen¬ 
tum distribution for the latter initial state, shown in 
Fig. S6(a), also exhibits a slight shift of the maxima to 
smaller absolute values of the peak position. One can 
understand this shift from the two-component picture 
discussed in Sec. IV of Ref. [63] as follows: We sep¬ 
arate the gas into the ballistically expanding unbound 
particles and inert doublons. The doublons form on time 
scales for which the cloud has not yet considerably ex¬ 
panded into an empty lattice [13]. Upon opening the 
trap, the doublons undergo the quantum distillation pro¬ 
cess [15, 64-66] that makes them accumulate in the cen¬ 
ter of the lattice, such that on the time scales probed 
by the experiment, they do not contribute to the qua¬ 
sicondensation (see the supplemental material of [13]). 
The number of doublons {d) = ANJ'^/[U{U — 6J)] can 
be estimated from second-order perturbation theory (see 
Eq. (39) in [63]). The condition for the position of the 
dynamically formed quasicondensate can then be formu¬ 
lated as E'/N' = —2J cos (q'd), where E' = — {d) U and 
N' = N — 2 (d) are the energy and the number of ballis¬ 
tically expanding particles, respectively. For U/J = 20, 
we find hq' = lA2{h/d)^ in very good agreement with the 
numerical data shown in Fig. S6(a). 

The main effect of a finite C// J = 20 is thus the slight 
shift of the peak positions in the dynamical quasiconden¬ 
sation of n{q). On the other hand, the shape of the peaks 
does not change considerably. Small residual potentials 
also cause a time-dependent shift of the position of the 
maxima [67], but do not lead to a broadening of these 
peaks. Our numerical results for HCBs (not shown here) 
are consistent with these predictions. 

From our previous work [13], we estimate that a small 
fraction of < 5% doublons is present in the initial state. 
Unfortunately, this regime of small doublon fractions can 
hardly be accessed with the tDMRG method as it would 
require large particle numbers, severely restricting the ac¬ 
cessible time scales. We have performed tDMRG simula¬ 
tions with one or two doublons in the initial state with an 
overall particle number of = 16. This translates into a 
fraction of doublons > 5% and we observe that in these 
cases (after averaging over many initial configurations) 
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the core expansion velocity Vc drops below the experi¬ 
mentally measured value, inconsistent with the results 
of Ref. [13]. It is nonetheless obvious from our tDMRG 
simulations shown in Fig. S6(b) that the presence of dou- 
blons in the initial state causes a significant reduction of 
the visibility of the finite-momentum peaks in n{q). Thus 
deviations between our experimental data and the results 
of the numerical simulations for HCBs can be partially 
attributed to doublons. 

Finally, we discuss the influence of the existence of 
some short-range coherence in the system present at 
tE = 0 on the expansion dynamics and the TOF distribu¬ 
tion. Such q = 0 short-range coherence is evident in the 
experimental data shown in Fig. 3(d) of the main text, 
and results in additional peaks in the TOF distribution 
beyond those expected from the emergent quasiconden¬ 
sates at finite momenta. This initial q = 0 coherence, 
which stems from an imperfect state preparation, clearly 
decays rather fast and one can see from two examples 
that it is not detrimental to the formation of quasicon¬ 
densates at finite momenta at sufficiently long expansion 
times. 

First, in the expansion from the correlated ground 
state of a Mott insulator at U/ J < oo, there is 
some short-range coherence present at tE = 0 already. 
Nonetheless, the results of Ref. [31], where the expansion 
from the correlated Mott-insultating states was studied, 
clearly show that sharp peaks at finite momenta form 
dynamically in n{q). Second, consider the sudden ex¬ 
pansion of hard-core bosons from a box trap, yet with a 
density below unity in the initial state and from a corre¬ 
lated state such as the ground state in that box. For this 
case, we show the time evolution of the quasimomentum 
distribution n{q) in Fig. S7. As expected, there is a max¬ 
imum in n{q) at g = 0 at tE = 0, which quickly dimin¬ 
ishes in favor of the quasicondensation peaks that emerge 
at hq = ±Tr/2{h/d). The quasicondensation peaks ulti¬ 
mately sit at the same position as if the expansion had 
started from the state given in Eq. (S14), but the distri¬ 
bution is not symmetric around hq = inj2(Hjd). Note 
that the reappearance of the peak at g = 0 observed at 
the longest expansion time t = 36t is a precursor to the 
dynamical fermionization occurring at very long expan¬ 
sion times [37, 38]. 


D. Current-carrying states and power-laws 

Both our expanding bosons and the melting domain- 
walls realize current-carrying states: In the first case, 
we are dealing with particle currents, in the latter case, 
with spin currents. Based on the scaling solution for the 
domain-wall melting [39], one can associate a current j (i) 
and a density n{x) to each point in space where x is the 
rescaled coordinate (see the main text). In the limit of 
large systems, each region described by the rescaled coor¬ 
dinate corresponds to an extended region in terms of the 
original spatial coordinate x, suggesting that one can use 



FIG. S7. Sudden expansion of HCBs from an initial state 
with an incommensurate density. We show n(g) for the ex¬ 
pansion from a box trap with A = 45 particles and an ini¬ 
tial density n = 0.9, plotted for expansion times t-E/r = 
0,7,14,22,36. 


the local density approximation to describe the system. 
In this picture, we seek homogeneous reference systems 
that have the same density n = n{x) and the same cur¬ 
rent j = j{x) as the expanding cloud at the rescaled 
position X. We then need two parameters to match these 
conditions. If we assume that the reference systems have 
periodic boundary conditions, then these two parameters 
are the chemical potential fixing the density, and the 
flux (j) through the ring, fixing the current. By numerical 
comparison we find that the chemical potential is linear 
in X while the flux (j) is independent of x and always given 
by 7r/2 (details will be published elsewhere). 

The question of power-law correlations can thus be re¬ 
formulated in two ways: First, do the (homogeneous) 
reference systems on which (j) enforces a nonzero current 
have power-law correlations and second, does this also 
apply to inhomogeneous systems. For simple ID spin sys¬ 
tems that allow a mapping onto free fermions and that 
have conserved currents, the answer to both questions is 
yes, according to the analysis of [39, 68]. The extension 
to interacting systems with or without exactly conserved 
currents or including integrability breaking terms will be 
addressed in a future publication. Examples for such 
systems are the spin-1/2 XXZ model or the ID Bose- 
Bubbard model. 

Beyond specific spin-1/2 chain models, related ques¬ 
tions pertaining to the existence of power-law correla¬ 
tions in excited states that carry a finite current have 
been addressed in various contexts, including highly ex¬ 
cited states in integrable ID systems [69, 70], quantum 
magnetism [41, 68, 71], hydrodynamics [72] and statis¬ 
tical physics [72-74]. Thus our work connects a broad 
range of branches of theoretical physics. 















